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Abstract. 

We study the quark-hadron phase transition within an effective model of QCD, 
and find that in a reasonable range of the main parameters of the model, bodies with 
quark content between 10 -2 and 10 solar masses can have been formed in the early 
universe. In addition, we show that a significant amount of entropy is released during 
the transition. This may imply the existence of a higher baryon number density than 
what is usually expected at temperatures above the QCD scale. The cosmological QCD 
transition may then provide a natural way for decreasing the high baryon asymmetry 
created by an Affleck-Dine like mechanism down to the value required by primordial 
nucleosynthesis . 
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1. Introduction 

It has been shown that the quark-hadron phase transition which occurred in the early 
Universe could lead to the formation of relic quark-gluon plasma objects, which survive 
today |L], 0]. Generally, it is admitted |1|, [3|, f|, |] that the transition occurred effectively 
at the critical temperature, which is of the order 100 MeV, the QCD energy scale. In 
that case, the quark content of the bodies which have been formed during the transition 
cannot be larger than 10 -8 M . Another possibility arises if the transition was delayed 
for some time. It then becomes possible that the quark plasma objects formed at the 
end of the transition appeared at a temperature much lower than T c , and more massive 
bodies may have been produced [§, [7|, j8|; these latter could then account for some 
fraction of the dark matter in the Universe. 

In this work, we perform a detailed analysis of the phase transition within an 
effective model of QCD, and we show that for the same critical temperature as what 
is usually thought, we can obtain a high degree of supercooling, which allows the 
possible formation of large "quark stars" with masses ranging from 10~ 2 to 10 M . 
We find that this is so because the Universe has grown exponentially during the quark- 
hadron transition, at a temperature T < 100 MeV, much lower than the temperature 
of classical inflationary models || . Moreover, we show that the exponential expansion 
is not balanced by so steep a drop of temperature, and thus increases significantly 
the total entropy of the Universe. This entropy production dilutes the density of any 
conserved or quasi-conserved quantity present before the transition, such as the baryon 
number. Therefore, our model requires a high value of these quantities at T > T c , as 
for example the large baryon asymmetry produced in some supersymmetric models [T(| . 



This leads us to incorporate the quark chemical potential in our model, and to study 
the cosmological quark-hadron transition in the context of both high temperature and 
quark number density. 

2. The Model 

Except in the quenched approximation, the theory of strong interactions QCD is 
too difficult, computationally, to give a reliable description of the phase transition 
between the high temperature and density phase of the quark gluon plasma and the low 
temperature and density phase of quarks bound in hadrons. In fact, lattice calculations 



1 1] do not seem able yet to really tackle the problem of finite quark densities. Because 
the light quarks have such a small mass compared with the transition temperature 
or chemical potential, we believe them to play an essential role in determining the 
quantitative features of the transition which is not taken into account in the quenched 
approximation, as shows for example the discrepancy between the values of the critical 
temperature computed with or without dynamical quarks. In considering the transition, 
we therefore work with a simple model in which non-perturbative QCD is replaced by 
an effective Lagrangian which incorporates the chiral symmetry of the theory. In this 
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model the quark fields are interacting with a chiral field formed with the n meson field 
and the scalar field a. The Lagrangian density is 

™ f - 11 

C = [Wk-fd^k ~ 9Mo- + it ■ n l5 )i(j k } + -d.ad^a + -d^n - V(a 2 + vr 2 ) (1) 



k=l 



which can be rewritten as |T2 

7lf 



1 1 



k=l 



where are the left- and right-handed components of the quark field ipk, U is an 



element of SU(2) defined by £U = a + ir ■ n and £ = (a 2 + 7r 2 ' 

The generalized self-interaction potential is the usual quartic function in £. We 
choose to write it in a seemingly complicated form, but such that the parameters f n , A, 
B are readily related to physical quantities. 

3" 



2 . 



1+3(4-] -4(4- 



/ \ fir 



(3) 



V(^) has its absolute minimum at £ = = 93 MeV, a value chosen to fit the observed 
pion decay rate. For this value, which corresponds to the physical vacuum, chiral 
symmetry is spontaneously broken; however, if n = isospin symmetry is preserved, 

A 2 f 4 

just as in nature. If B < — jyS the potential has a second local minimum at £ = 0, 
corresponding to a chirally symmetric, metastable "false vacuum" with the energy 
density B. This is analogous to the perturbative vacuum of the MIT bag model [13], 
and B can be interpreted as the bag constant in that model. Various versions of the 
previous model have often been used to study the low energy hadron spectroscopy (T^| or 
heavy ion collisions [I5|, and phenomenological fits to light hadron properties give 

#1/4 

between 100 and 200 MeV f[6| . Small oscillations of £ about the minimum at £ = f w 
correspond to a scalar particle of mass = Xf n and three massless pseudoscalar pions. 
We can notice that the £ field is a chiral singlet, and the associated scalar particle may 
be interpreted as representing the condensate arising from non-linear interactions of the 
gluons (glueballs). As the lightest glueballs are believed to have a mass in the range 
1.5 to 1.7 GeV JT7L P|, A is between 16 and 19. The coupling constant g gives mass to 
the quarks, and thus helps to break the chiral symmetry; following the fits to hadronic 
physics, we take it to be larger than 10, so that quarks have an effective mass larger 
than 1 GeV in the physical vacuum: it is then energetically unfavourable for them to 
exist in the phase with £ = f n . However, the actual value of g is not important for the 
following calculations. 

In order to study the transition which took place in the early Universe, at a 
temperature of the order 100 MeV, we must implement finite temperature in our model. 
According to the common lore, it seems irrelevant to take into account a possible 
chemical potential, because the quark density in the cosmological context is too small 
to be significant (see e.g. |I9|, pp. 530-531). To this effect, we must add to the potential 
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given by equation (P) contributions like [^(], |21 



ff(T)= -2T| I 01n(l + e-^), (4) 

E(k) = ^k 2 + g 2 e 
for each fermionic degree of freedom (spin, flavour, colour), and 

A(r)=T /(S 5in ( i - e " E ' ( " /T )- (5) 



E'(k) = Jk 2 + 



2 



for bosonic degrees of freedom. These terms come from the one loop approximation and 
represent the effect of the thermal excitations of the quark-antiquark pairs and of the £ 

/„ 2 „\ 1/2 

field. These latter have a mass ( ^p- j larger than in both phases with £ = and 

£ = f n , and are strongly suppressed by the Boltzmann factor e~ E ' l T at T < 100 MeV. 
The same arguments afford us to discard the excitations of qq pairs in the physical 
vacuum £ = /„., but not in the false vacuum £ = 0, where the contributions given by 
equation (§J) lower the value of the effective potential. For £ close to 0, which implies 
i<l, the integral can be expanded as |20, |22|] 



T 

f ( T ) = _!![! T 4 + 1^2^2 T 2 (6) 

Jn 1 360 24 y ^ y 1 

Actually, we have checked numerically that equation @ represents a very accurate 
approximation. We include in our model the three light quarks. Each contribute for 3 
degrees of freedom of colour, and 2 of spin, giving 18 degrees of freedom in all. The u 
and d quarks are considered massless, whereas for the strange quark we will keep the 
quadratic term of equation (|6]) to take into account its finite mass. 

With these modifications, the self-interaction potential at finite temperature and 
£ = can be rewritten 

V T (£) = B- a T T 4 + 7T T 2 (7) 

where a T = ^ and j T = \m 2 s1 m s = (60 - 170) MeV.pl 



20 I 1 4 

A sketch of Vr(£) for different values of the temperature is given in figure [I]. It has 
two minima, at £ ~ f n with the same value as at zero temperature, and at £ = 0, where 
the value of the potential is lowered. The overpressure, difference in free energy density 
between these minima, is 

AP = B + lT T 2 - a T T A (8) 

This difference vanishes when 



t =t c = \/ 7t+v 2 7 ; t +4qtjB 0) 

For this value, both minima of Vr(£) are degenerate, which gives a first order phase 
transition. Even though the order of the QCD transition is still a highly debated issue, 



Possible Cosmological Implications of the Quark-Hadron Phase Transition 



5 



this is in agreement with previous theoretical predictions using the linear a model [3J 
which, incidentally, is a particular case of our model when A 2 / 4 = 8B, or with recent 
indications from lattice calculations with 3 degenerate Wilson quarks or with 2 massless 
quarks and a light s quark fS3 |. T c is then the critical temperature of the transition. If 



T > T c , the minimum at £ = becomes the absolute minimum of Vr(£). 

Besides the quarks, we first include in our model a thermal bath of relativistic 
particles which influence the transition only through gravitational effects: at this stage, 
they represent the "physical vacuum". The contribution to the free energy density 
of both phases of these spectator particles (photons, electrons, muons and neutrinos) 
is given by expressions (|J) and (j|) for fermions with their antiparticles and bosons 
respectively. These terms can be expanded as ff = — ||^-T 4 + ^m 2 T 2 , with a mass term 

only for the muon, and /& = — fjjP 4 for the massless photon, so we get an extra 
contribution 

/« = -^f^ 4 + T2< T2 = ^ + (10) 
contributing to both phases. However, if we consider the overpressure AP between the 
phases with £ = and £ = f w , this extra free energy density cancels, and AP is still 
given by equation (|8|). 

In a more exhaustive calculation, we will also include in the hadronic phase £ = f n 
the lightest hadrons, the pions, and show that their influence on the results, which we 
discuss later, is not qualitatively important. 

3. Bubble nucleation. Dynamics of the Universe 

After the Big Bang, the temperature of the early Universe is higher than the critical 
temperature of the quark-hadron transition (equation ^): the phase with £ = is 
the more stable. The massless quarks are deconfined - this is the so-called quark- 
gluon plasma -, and chiral symmetry is preserved. According to the standard model of 
cosmology, the Universe expands, while its temperature decreases. When it drops below 
T c , bubbles of the true vacuum £ = /„., which is now the stable phase, begin to appear 
within the quark plasma. If we consider only homogeneous nucleation, the nucleation 
rate per unit volume is given by |25 



r(T) = CT 4 (^ 3/2 e-^ (11) 



C is a multiplicative coefficient of the order unity jSJ|. For a spherical bubble, S3 is 
the stationary value of the functional F = 4tt J r 2 [|(^) 2 + Vr(0] dr. In the thin-wall 
approximation [p5|| , this expression can be replaced by F = — |7rrgAP + 4ivr 2 s, where 
r is the bubble radius, the overpressure AP is the difference in Vr{£) between the two 
minima, given by expression (H), and s is the surface tension. This latter is 

'"VW&)d£ = \fil{y) (12) 
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with 

I(y) = f (l-w)y^ 2 +|(l + 2u-3<) du, V = ^i- (13) 

The graph of I(y) in figure shows that | < I(y) < |. It should be noted that s 
is entirely fixed by the parameters of the zero-temperature model, determined at low 
energy. Its value is then 50 < s < 120 MeV/fm 2 , higher than the estimates given by 
lattice computations without dynamical quarks p7j , but in agreement with the surface 
tension obtained by Burakovsky pS] within an effective model of QCD different from 
ours. 

F is stationary at a radius r = and takes the value 

= 16vr^_ 16tt [Xmy)} 3 
3 3 AP 2 3 AF 2 1 ' 

Let us call t c the time at which the temperature equals T c . When t > t c , both 
phases with £ = and £ = f n coexist. If is the fraction of the Universe occupied by 
the physical vacuum, and therefore 1 — x(t) the fraction taken up by the quark plasma, 
the energy density of the matter which fills the Universe is 

e(t) = [l-x(t)]e i=0 (t)+x(t)e^ f7v (t) 

= [1 - x(t)} [e q (t) + e v (t)\ + x{t)e v (t) (15) 

e q , e v , e^=o and £5=/^ stand respectively for the energy densities of the quark-gluon 
plasma, the thermal bath of photons and leptons, the metastable vacuum and the 
physical vacuum. They are given by 

e q = B- lT T 2 + 3a T T 4 (16) 
e v = 3a v T* - lv T 2 (17) 

£§=0 = e q + £ v 

Expression (|l^) can be rewritten 

e(t) = [1 - x(t)]s q (t) + e v (t) (18) 

It appears as the sum of two contributions corresponding to the quark plasma, which is 
going to disappear during the phase transition, and to the background of non-interacting 
particles. If we replace the energy densities e q {t) and e v {t) with their temperature- 
dependent expressions (|TE| , |TTD , we obtain 

e{t) = [1 - x(t)} [B - j T T(t) 2 + 3a T T(t) 4 } + 3a v T(t) 4 - lv T{t) 2 (19) 

If we now suppose that the nucleation of the bubbles of phase with £ = f n is 
isotropic enough, the Universe is still described by the Robertson- Walker metric. We 
can then deduce from Friedmann's equation and the previous expression the equation 
which governs the expansion : 



1 dR , . / 8tt . . . . 

(*) = W^-2- £ (*) ( 20fl ) 



R(t) dt y ' V 3m 2 P 
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dR 
~Rdt 



in 



3mp 



[(1 -x){B- 7T T 2 + 3a r T 4 ) + 3«„T 4 - lv T 2 } 



(206) 



where mp is the Planck mass while x and T are time-dependent. 

In the beginning, few bubbles appear, and the temperature still decreases. As time 
goes on, the size of bubbles increases, because of both the propagation of the bubble 
walls and the expansion of the Universe. For example, if we consider a bubble of radius 
r(t) which appeared at time £j, its growth is given by 

" dt' 



dr 
dt 



v + ?*4^|, hence 



r(t) = vR(t) 



R{t>) 



(21) 



v is the propagation speed of the walls, which we shall take equal to the velocity of light 
in our calculations. 

While the size of existing bubbles increases, new bubbles form. Thus the number 
density of nucleation sites at time t is 



N(t) 



[1 - x(i/)} T(t>) 



R{t) 



df 



(22) 



This expression shows that the number N[t) [-R(t)] 3 of bubbles in a unit comoving volume 
which exist at t, i.e. those which appeared at any time t' < t, is increasing. However, 
the number of bubbles per unit physical volume N(t) may either increase, when the 



R(t>) 



1 3 



becomes 



nucleation rate T(t') is large, or decrease, when the expansion factor 
large. In other words, the mean physical distance between nucleation sites may either 
decrease, in the first case, or increase as the Universe expands. Nonetheless, the fraction 
x(t) of the Universe filled up with physical vacuum is growing with the number of 
bubbles. As we will see later, it is necessary to take into account the expansion of the 
Universe in the evolution of this fraction, so that 



x(t) 



e W (-j\l-x(t')}V(t')^ vR(t)J 



dt" 

W) 



dt' 



(23) 



This expression is analogous to the formula derived by Guth and Weinberg |2!| in the 
context of GUT phase transitions, with the exception of the extra factor 1 — x(t'), which 
takes into account the fact that hadronic bubbles can only appear within the quark 
plasma. 

The evolution of T during the phase transition is easily determined by using the 
energy conservation, which reads 







s(t) is the energy density, given by expression (0), whereas the pressure is 
P = (1 - x)(a T T 4 - 7r T 2 -B) + a v T 4 - lv T 2 



(24) 



(25) 
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The time-dependences of e or P are connected with those of T and x, so we get 
dT _ (1 - x) {6a T T 3 - 3j T T) + 6a v T 3 - 3^ V T dR 
dt (1 — x) (6arT 2 — jt) + 6ct v T 2 — j v Rdt 

3a T T 4 - 7T T 2 + B dx 

(1 - x) (12a T T 3 - 2 7T T) + \2a v T 3 - 2 lv T dt [ ' 

^ is the sum of two terms : the first one is the usual contribution of the expansion 
of the Universe, which decreases the temperature. The second contribution comes from 
the replacement of the false vacuum £ = with the physical vacuum £ = f n , which has 
a much lower specific heat, and corresponds to the release of latent heat, which tends to 
increase T. In earlier studies fl], |3|, §|, §J, it is assumed that the latent heat released by 
the first bubbles of hadronic phase formed when the temperature drops under T c quickly 
reheats the Universe up to the critical temperature, and further bubble nucleation is 
suppressed. For the remainder of the transition, i.e. some microseconds, the temperature 
is kept constant by the latent heat due to the growth of the hadronic bubbles, which 
balances the expansion of the Universe. This latter remains negligible: equation 
reduces to r(t) m v(t — t c ), while the expansion factor in expression (^) is almost 
equal to 1. In this work, we investigate this point more carefully, following the detailed 
evolution of the temperature. In particular, we will allow out-of-equilibrium coexistence 
at a temperature different from T c between the quark plasma and the hadronic phase. 
For example we can expect from equation (|26|) that the product R(t)T(t) will not 
remain constant during the transition, and the total entropy of the Universe, which 
is proportional to (RT) 3 , will not be conserved. 

In many calculations concerning first order cosmological phase transitions |29], SDj , it 
is assumed that R(t)T(t) = constant during the transition, i.e. that entropy is conserved. 
In that case, the second term of equation (|26l) is neglected, and following the work of 
Coleman |3T| it is thought that all the latent heat is taken by the bubble walls. At the 
completion of the transition the Universe would then reheat as the energy in the walls 
is dissipated by their collisions fl32"| , and T would quickly increase up to a temperature 
close to its initial value T c . In this work we will see that there is no such steep drop 
of temperature, but that the latent heat is released gradually, and holds T rather high, 
though not necessarily at T c , during the entire transition. When this latter is completed, 
the temperature may then be significantly smaller than the critical temperature T c . 

After a while, bubbles begin to coalesce and to enclose lumps of the phase with 
£ = 0. The quark number enclosed in such a chunk is conserved, whereas the size of the 
lump decreases as the bubble walls propagate, and therefore its density increases, until 
the pressure of the quark plasma trapped inside balances the pressure of the physical 
vacuum £ = f n . At that time, we can consider that x ~ 1, because the regions filled 
with quark plasma are but a minute fraction of the Universe, and the transition has 
come to an end. Let us call tf the time at which x(tf) = 1, and N(tf) the corresponding 
number density of nucleation sites. To simplify the calculation we have not considered 
the dispersion in the number of trapped quarks in the chunks as given by the statistical 
distribution of nucleation sites. To estimate the mean number of quarks per body, we 
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shall assume that these sites are on the vertices of a cubic lattice. At the centre of each 
cube is a quark clump. If no phase transition had occurred, the quarks enclosed would 

r 1 3 

take up a volume V} = N * t y which corresponds to a physical volume V c = Vf j^j^ 
at the beginning of the transition. Therefore, the body contains M q = n q (t c )V c quarks, 
where n q (t) is the quark number density at time t, which can be easily related to 
the photon number density n 7 = ^C(3)T 3 and the quark to photon ratio — . The 
conservation of the quark number when there is no baryon number violating interaction 
reads n q (t c )[R(t c )] 3 = n q (tf)[R(tf)} 3 , so it is possible to rewrite Af q = n q (tf)Vf = -jfj^, 
hence 

^=-tC(3)^(*/)tSt (27) 

— (tf) is the quark to photon ratio at the end of the transition. We suppose that 
this value remained unchanged till the beginning of the Big Bang Nucleosynthesis at 
T ~ 1 MeV, and therefore take in our calculations ^ = 3 -^ ~ 10~ 9 , with ng the baryon 
number density As a matter of fact, this value of ^ is valid in the hadronic phase, 
not in the quark plasma. However, estimates JI], ^| show that the baryon number density 
in the quark phase should be larger than in the hadron phase, and thus our assumption 
tends to underestimate Af q . 

4. Results 

We have tested our model with different values of the parameters B, m s and A. With a 
high value of B, such as B l / A = 200 MeV, and m s = 100 MeV, the critical temperature 
is T c ~ 148 MeV. If we take A = 18, we find the usual results [f[], EL although our model 
is different, and the number of quarks trapped in a nugget at the end of the transition, 
which is completed in some microseconds, is M q — 10 40 . In that case, the standard 
model of cosmology shows that the horizon at tf is too small to allow any object with 



a mass exceeding 10 M Q to be formed [22, 34], which is in agreement with what we 



find. However, it may be worth noting that the upper limit given by this simple but 
powerful argument will no longer be valid if the Hubble volume is larger than what is 
given by the strictest standard model of cosmology (i.e. R(t) oc t 1 / 2 ). 

Another limiting situation appears when B takes a very low value, as for example 
5 1 / 4 = 50 MeV. In this case, we find that x(t) cannot reach 1, and therefore the 
transition never ends. The energy density which confines quarks within solitons is so 
weak that few bubbles of the physical vacuum £ = f n are nucleated, and they are so 
far away from each other that their walls cannot meet while the Universe expands. An 
analogous phenomenon occurs when the surface tension s is too high. This happen also 
if the mean distance between bubbles [-V(tj)] 1 / 3 is smaller i.e. if the value of B is larger, 
but if simultaneously the speed of the bubble walls v is decreased. Once again, the 
expansion of the bubbles cannot catch up with the expansion of the Universe. 

Now if B takes intermediate values, like B 1 / 4 = 120 MeV (i.e., in other units, 
B = 0.61 fm~ 4 , which is obtained from fits to hadronic properties in [35fl) with 
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m s = 100 MeV, which gives a critical temperature T c = 90.1 MeV, and A ~ 17.77, 
corresponding to a glueball mass of 1.6 GeV, the number of quarks enclosed is M q — 10 57 , 
i.e. a quark content of approximately 0.3 M . We have plotted the evolution during the 
transition of several quantities in such a case: the fraction x(t) of Universe filled with 
physical vacuum, the temperature T(t) and the scale factor R(t). 

We can see in figure [5] that x(t) remains close to for some time, when few bubbles 
have yet appeared. Then it rises quickly, though not instantaneously as in the "fast" 
scenario, and tends towards its asymptotical value 1. The time scale of the transition 
is now of the order of the millisecond, rather than the microsecond. 

The temperature (figure |4j) behaves as we had foretold. In the beginning, it 
decreases with a law close to the t~ 1 / 2 law of the standard model: the first term of 
the r.h.s. of equation (|26"D is dominant. Then it increases, though not much, when one 
phase is replaced with another: the influence of the term in which represents the 
release of latent heat, is the more important; we can check that the fast growth of x 
and the increase in T are simultaneous. Finally, T decreases like t~ 1 / 2 again, when the 
Universe is mostly filled up with physical vacuum, down to Tf = 17.7 MeV. 

The scale factor R is plotted in figure || Obviously, it does not follow the t 1//2 law 
of the standard model, but rather seems to grow exponentially, to end up with a value 
Rf ~ 7 x 10 3 i? c . It should be noted that in this case, the dilution factor in equation 



R(t) 



l 3 

can be as small as 10 -11 ! During the transition, the 



is no longer negligible, since 
scale factor has been multiplied by ~ 10 4 , and therefore the Hubble radius is 10 4 times 
the radius predicted in the strict standard model: the Hubble volume is ~ 10 11 times 
larger than what is expected, and there is no contradiction between the number M q of 
quarks enclosed in a nugget and the total number of quarks within the Hubble volume. 

If we compare our model with Guth's inflationary model ||30|| , we can trace whence 
such a growth comes. As in his case, when the temperature has dropped sufficiently, its 
contribution in our equations (eq. [19]) can be neglected with regard to the metastable 
vacuum energy density B. This latter then drives the expansion However, in 

our model the exponential growth comes to a natural end according to equation ( J206p , 
since the factor 1 — x(t) vanishes when the transition is completed. A second difference 
between our model and old inflation a la Guth regards entropy. In his model entropy 
is conserved during the exponential expansion and increases, due to reheating when 
bubbles collide, only at the completion of the transition. On the contrary, here entropy is 
constantly increasing during the quark-hadron phase transition, as the product R(t)T(t) 
does not remain constant. To be more precise, the scale factor R is multiplied by 7 x 10 3 , 
whereas T only decreases from 90 to 18 MeV: the product RT is multiplied by more 
than 1.4 x 10 3 , and hence the entropy increases by a factor 3 x 10 9 . 

This increase in entropy has radical consequences on the evolution of the ratio — , 
which is a measure of the baryon asymmetry of the Universe, during the transition. We 
have already seen that if there is no baryon number violating interaction, the quark 
number density n q decreases as R~ 3 , and that n 7 oc T 3 ; hence, the ratio — oc (RT)~ 3 is 
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proportional to the inverse of the entropy. In particular, if, as we have found before, the 
entropy is multiplied by 10 9 during the quark-hadron transition, the baryon asymmetry 
is divided by the same factor. Therefore, if at the end of the transition — ~ 10 -9 as 
required by primordial nucleosynthesis, it means that at T > T c the ratio was of the order 
1. The results found here contrast sharply with normal adiabatic expansion in which 
this baryon asymmetry would not change. However, when ^ ~ 1, the approximation 
made so far of neglecting the quark density or the chemical potential of the quarks with 
regard to T is no longer valid. 



Uf(T, (i) = -T 



(2f 



5. Finite chemical potential 

To handle the possibility of large quark densities, we must now incorporate chemical 
potential into our model. Within the one loop approximation, the contributions of 
fermionic degrees of freedom which should be added to the effective potential V(£) at 
zero temperature and zero chemical potential are no longer given by expression (^), but 
by 

7 (0 ln {1 + e- <£, "-' ,/T ) + / (gs >« (i + ^> E ^n 

The first term is the fermion contribution, the second is from the antifermions. We 
have changed the notation from ff to Uf because at finite /i, the free energy and the 
thermodynamic potential do not coincide any more. For the same reason as before, the 
excitations of the scalar field vanish at both minima of V(£), whereas those of quark- 
antiquark pairs are important only for £ close to 0, where it becomes possible to expand 
the right-hand side of equation ( p8|) as 

ujf(T, n) = -— T 4 - -fT 2 - -L/i 4 + -g 2 eT 2 + A^CV (29) 
a >w 360 12 p 24?r2 p 24 y s g7r2 y s f* v ) 

As previously, we keep the mass terms only for the strange quark. The self-interaction 
potential at £ = reads 

VtAO = B - a TT 4 - I^t^T 2 - + lT T 2 + 1 ^ 2 (30) 
7vr 2 3 3 1 2 3 2 

As before, we consider B and m s as parameters, while \x is related to the quark number 
density by n q = 4a^ 3 + 2-f^fiT 2 - 

Therefore, the difference in pressure between the two minima of Vr )M (£) is 

AP = B + lT T 2 + 7 M /i 2 - a T T 4 - lilT ^i 2 T 2 - (31) 

This in turn changes the expression of the critical temperature, given by the condition 
AP = 0, and of the bubble nucleation rate (|TTD. 

The energy density which drives the expansion of the Universe through Friedmann's 
equation ( |20 oj) is now 

e = (1 - x) (B - -f T T 2 - 7/ y + 3a T T 4 + 3 lflT ^ 2 T 2 + 3^/i 4 ) + 3a,T 4 - lv T 2 (32) 
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whereas expression (j26|) , which gives the evolution of temperature during the transition, 
becomes 



dT 3T 



^i-x) ( 2ar r 2 + w 2 - 1T - ^ggg^ ) + 2a,r 2 - lv 

B + 3a T T 4 + 3c^ 4 + ^^tfT 2 - j T T 2 - dx 
1 - x) (Qa T T 2 + lflT ^ - 7r - 6 ^y;^- 7 J + 6 «- T2 - "ivT 2 



(33) 



The remainder of the calculation of the number M q of quarks trapped in a typical 
chunk is changed accordingly. However, whereas the number of the model parameters 
B, m s and A remains the same as in the vanishing chemical potential case, there is 
now an important extra requirement, namely that the quark to photon ratio after the 
transition should be equal to 10~ 9 , in agreement with the value required for primordial 
nucleosynthesis. From now on, we choose to take 5 1/4 = 120 MeV and m s = 100 MeV, 
whereas A will be determined by the above condition. In other words, once 2 of our 3 
parameters have been fixed, for example the energy density of the false vacuum and the 
mass of the strange quark, the value of the third parameter cannot be chosen at will, 
but is dictated by the amount of entropy which has to be released during the transition. 
Hopefully, as we will see, the "good" values of A come exactly in the range which gives 
the glueball mass. 

If we take a chemical potential ^(t c ) = 31.44 MeV, equivalent to — (t c ) = 1 at 
the critical temperature T c = 89.96 MeV, and with A ~ 17.77, at the end of the phase 
transition, the number of quarks enclosed in a chunk is M q — 2 x 10 57 , i.e. approximately 
0.5 M . The final temperature is Tf = 17.91 MeV and the ratio of the final and initial 
scale factors ~ 4 x 10 3 , so that the quark to photon ratio after the transition has 
become ^ ~ 2x 10~ 9 . However, we have to wonder how an important quark number 
density could have arisen in the early Universe, at temperatures above the QCD scale. 



Affleck and Dine |10[ have proposed a mechanism for baryogenesis which could give 
such a high initial value of — . In their model, the decay when supersymmetry breaking 
effects become important of the large vacuum expectation value of a scalar fermion, 
acquired either by quantum fluctuations at the Planck epoch or after an episode of 
inflation [37], leads at a temperature of the order 10 TeV to a baryon asymmetry which 



can be as large as — = 10 3 . The actual value of the asymmetry depends on several 
parameters of the model, such as the expectation value of the squark field which decays 
or the CP-violating phase required to produce more quarks than antiquarks. After its 
production, this huge baryon asymmetry could still have been wiped out, long after the 
SUSY breaking epoch, by the so-called sphaleron transitions |38[] arising, for example, 



from the electroweak anomaly. However, it was shown in reference |55 that a baryon- 
lepton number (B — L)[| excess can be produced by the Affleck-Dine mechanism along 
with the usual B excess. This is sufficient to prevent any electroweak sphaleron induced 

f Throughout this paragraph, B denotes the baryon number, not the energy density of the chirally 
symmetric vacuum £ = 
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Table 1. Values of Af q with B 1 / 4 = 120 MeV, m s = 100 MeV, A variable. 





A 


T 

± c 

(MeV) 


(MeV) 


(MeV) 


*/ 

(ms) 


Rf 
Rc 


£(*/) 






0.01 


17.77 


90.10 


8.98 


19.81 


4.612 


727 


2 x 10" 


-9 


1.8 x 10 55 


0.1 


17.77 


90.10 


15.48 


18.56 


5.990 


1739 


2 x 10" 


-9 


1.8 x 10 56 


1 


17.77 


89.96 


31.44 


17.91 


7.391 


3974 


2 x 10" 


-9 


1.8 x 10 57 


10 


17.56 


87.28 


64.94 


17.60 


8.706 


8405 


2 x 1(T 


-9 


1.5 x 10 58 


20 


17.27 


83.70 


78.41 


17.64 


8.636 


10 4 


2 x 10" 


-9 


5.2 x 10 57 


30 


16.99 


80.25 


86.05 


17.83 


7.970 


10 4 


2 x 10" 


-9 


1.6 x 10 57 


40 


16.73 


77.09 


90.98 


18.78 


6.190 


5523 


2 x 10" 


-8 


3.1 x 10 56 



B violation which would erase the above mentioned baryon asymmetry. Furthermore, if 
we leave aside the electroweak symmetry breaking transition, which, within the standard 
model, does not seem to affect much the baryon asymmetry | iD| , and assuming normal 
adiabatic, radiation dominated evolution, the quark to photon ratio left by the Affleck- 
Dine baryogenesis should not be significantly changed down to the QCD transition 
considered here. 

Therefore, we have taken several different ratios — (t c ), keeping in mind the 
necessary requirement ^(tf) — 10~ 9 at the end of the transition. This requirement 
can be met, for moderate values of the baryonasymmetry, by changing slightly the value 
of A. However, for larger values of ^(t c ) the change in A is more pronounced, in order 
to get enough expansion to dilute the initial baryon asymmetry to the desired value. 
The results with B 1 ^ = 120 MeV, m s = 100 MeV and different values of A are shown 
in table |l|. 

The first result which is striking in this table is the apparent proportionality between 
the initial quark number asymmetry — (t c ) and the number of quarks enclosed in a 
typical quark star. When — (t c ) grows from 0.01 to 10, Af q takes values corresponding to 
a quark content of 0.005 to 5 M . This seems natural: let us assume that the parameters 
B and m s (remember that A is no longer a free parameter) determine the geometrical 
repartition of bubbles, i.e. the mean distance between nucleation sites. Then, the more 
quarks you put initially in the volume which corresponds to a future quark object, the 
more quarks you will have in that object in the end. However, when — (£ c ) is larger 
than 20, the number of quarks M q does not increase anymore, but on the contrary it 
decreases. This means that it is impossible to form quark stars with gigantic masses, 
even with fine tuning, whereas in previous studies such fermion soliton stars could have 



a mass as high as 10 M [[16], [y]] 



So far, we have assumed that the quark-hadron phase transition takes place between 
the quark phase and the QCD vacuum. In a subsequent study, we have considered the 
more realistic case with the transition between the quark plasma and a hadron phase 
containing the lightest hadrons, the pions. These contribute only to the phase £ = f n , 
by an additional term — |qT 4 + ^-T 2 to the self-interaction potential Vr, M (£ = fw), 
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which changes the expression of the overpressure (eq. pTl) and the value of the critical 
temperature. Similarly, a term x(j^T 4 — ^-T 2 ) must be added to the energy density 
(eq. [32]), thereby modifying Friedmann's equation and the energy conservation (eq. |33|) . 
The resulting calculations show that the qualitative features presented above remain 
unchanged, provided we slightly tune the value of A. For example, with _B 1//4 = 120 MeV, 
m s = 100 MeV and an initial baryon asymmetry — (t c ) = 10, the baryon asymmetry at 
the end of the transition equals 10 -9 if we take A = 17.76, and the number of quarks 
enclosed in a typical quark star is Af q = 1.7 x 10 58 , close to the value obtained when 
pions are neglected. 



6. Conclusion 



We have studied the quark-hadron phase transition within an effective model of QCD, 
the chiral quark model, including finite temperature and chemical potential. We would 
like to stress that this type of effective models is widely used in various fields of physics 
like the low energy hadron physics or relativistic heavy ion collisions. It is also worth 
noting that the so-called inflaton or dilaton scalar field models often used in cosmology 
are special cases of this model. In our work, all the parameters of the model are fixed 
by fits to the static properties of hadrons, and we have seen that in a reasonable range 
of the key parameters B, m s and A, quark plasma objects with very different masses 
can have been formed at the end of the transition, which are likely to survive until 
our epoch if their quark content is M q > 10 40 (see the discussion and the references in 
|34f). The most interesting case regards the formation of a body with stellar mass, at a 
temperature much lower than the critical temperature of the transition. Such quark stars 
could be identified with some of the dark bodies detected by microlensing experiments 
in the halo of our galaxy [fE^| . A new feature of our calculations, independent of the 
size of the formed quark bodies, is the presence of a period of exponential expansion 
of the Universe during the transition which ends in a natural way with the transition. 
This growth, unexpected at such a temperature, could account for the observed baryon 
asymmetry in the Universe, if a high value of the quark to photon ratio had been left 
before the transition, at T > 100 MeV by a baryogenesis mechanism d la Affleck-Dine, 
i.e. if we drop the usual assumption that fermion densities in the (very) early Universe 
were negligible: the entropy production by the released latent heat could have diluted 
the asymmetry down to the value it must have at T < 10 MeV to agree with the 
standard Big Bang nucleosynthesis. Actually, various models of entropy production 
have been designed for this purpose at very high energies by several authors flclfl , but 
our dilution mechanism is different from those earlier proposals in that it involves only 
the physics at the QCD scale. 
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Figure 3. Evolution of the fraction of the Universe filled up with the phase £ 
during the transition in the case B = (120 MeV) 4 , m s = 100 MeV, A ~ 17.77. 
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Figure 4. Evolution of temperature during the transition. 




Figure 5. Evolution of the scale factor during the transition. 



